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ABSTRACT 


The  report  is  concerned  with  the  resonant  behaviour  of  general 
structural  systems  under  single  frequency  harmonic  excitation.  In 
particular  it  is  concerned  with  the  parametric  and  nonlinear 
phenomena  which  arise  when  there  are  certain  integral  relationships 
between  the  excitation  frequency  and  some  of  the  natural  frequencies 
of  the  structure  (external  resonance) , and  amongst  the  natural 
frequencies  themselves  (internal  resonance). 

Qualitative  experimental  evidence  is  presented  illustrating  the  I 

i 

occurrence  of  this  type  of  behaviour  in  simple  structures  having  up  to 
four  distinct  modes. 

In  the  theory  the  equations  of  motion  retaining  linear  parametric 
terms  and  including  quadratic  inertial  nonlinearities  are  considered. 

A general  theoretical  approach  is  presented  which  is  applicable  to  any 

i 
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multimode  structure.  The  theory  predicts  the  resonance  conditions,  j 

and  in  the  vicinity  of  each  yields  simpler  equations  giving  an  c 

approximate  solution.  Generally  these  equations  require  a numerical  i 

solution  for  a given  structure.  In  certain  cases  further  analysis  on 
a reduced  form  of  the  equations  yields  analytical  expressions  which 
define  the  boundaries  of  the  resonance  region  in  the  excitation 
amplitude- frequency  plane. 

The  theory  is  applied  to  lumped  mass  mathematical  models  (based 
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PRINCIPAL  NOTATION 


Aj,(t) 


Excitation  amplitude  (displacement) 

Slowly  varying  modal  amplitude 

D.C.  (zero  frequency)  slowly  varying  amplitude 


Inertia  matrix 


FN1,2 


Stiffness  matrix 


= ” ^f^  ’ linear  forced  response  component 

Excitation  amplitude  (force,  normal  coordinates) 
Excitation  frequency  (computer  diagrams) 

Response  frequencies  (computer  diagrams) 

Natural  frequencies  (computer  diagrams) 

Beam  length 

Nonlinear  coefficients  (quadratic) 


SCI,  2 


Nonlinear  coefficients  (quadratic) 

= f(5,^^j,  Up j j > ^j)'  Autoparametric  terms  in  variational  equations 
Normal  coordinate 

= f(K^.  ,£^..  ,m^.j^,fij^)  parametric  type  terms  in  the  variational 
^ equations 

Modal  matrix 

Amplitude  scaling  factors  (computer  diagrams) 


Negative  damping  coefficient 


= (fl^  + - (jjp  , detuning  parameter 


= • A.  D.C.  additions  to  variational  equations 

nr  lo  r’  ^ 

1 = 1 


Small  parameter 


Parametric  term 


Beam  stiffness 


Slowly  varying  modal  phase 
Forced  response  phase 
Natural  (linear)  frequency 
Response  frequency 
Excitation  frequency 
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INTRODUCTION 
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1.1  GENERAL 

This  report  is  concerned  with  some  aspects  of  the  'resonant' 
response  of  structures  in  situations  in  which  parametric  and 
nonlinear  influences,  which  are  always  present,  become  of  considerable 
significance.  In  'uch  situations  response  estimations  using  the 
usual  constant  coefficient  linear  modelling  of  the  structure  are 
quite  useless.  The  'smallness'  of  the  nonlinear  terms  or  the  time 
varying  coefficients  is  not  sufficient  to  prevent  them  from  having 
an  overwhelming  effect  on  the  response  in  the  course  of  time.  In 
many  engineering  situations  the  excitation  is  applied  to  the  structure 
In  a more  or  less  steady  oscillatory  fashion  for  an  appreciable 
duration  thus  allowing  the  development  of  the  influence  of  these 
effects. 

In  more  precise  mathematical  modelling  of  structures  in  dynamic 
conditions,  nonlinear  quadratic  and  higher  order  inertial  terms  are 
nearly  always  present.  The  adjective  'inertial'  is  applied  to  these 
terms  which  have  their  origin  in  the  kinetic  energy  T of  the  structure. 

The  existence  of  quadratic  nonlinearities  from  this  source  is  very 
important  but  curiously  it  has  received  little  attention  in  the 
literature.  Elastic  nonlinearities  have  received  much  more  consideration 
but  they  are  usually  a symmetric  form  so  that  they  appear  as  odd 
functions  through  cubic  and  quintic  nonlinear  terms. 

This  report  concentrates  on  behaviour  characteristics  which  arise 
from  or  can  be  explained  by  the  existence  of  the  quadratic  inertial 
nonlinearities.  The  higher  orders  of  nonlinearity  produce  their  own 


phenomena  and  as,  in  certain  excitation  frequency  bands,  the  linear 
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model  has  to  be  augmented  by  consideration  of  quadratic  nonlinearities 
so  this  model  in  turn  may  also  require  cubic  nonlinearities  to  be 
included  to  produce  qualitatively  certain  phenomena  or  quantitatively 
a required  degree  of  accuracy.  Fortunately  it  seems  to  be  the  case 
that  in  general  the  influence  of  higher  and  higher  orders  of 
nonlinearity  of,  presumably,  smaller  and  smaller  magnitude  becomes 
virtually  undetectable  in  the  presence  of  damping  so  that  an 
examination  of  terms  of  higher  order  than  the  cubics  is  seldom 
necessary  to  explain  observed  behaviour. 

In  structural  dynamic  analysis  'ordinary'  forced  excitation,  an 


independent  time  function  constituting  a nonhomogeneous  'right  hand 
side'  to  the  equations  of  motion,  is  usually  the  only  forcing  action 
considered.  However  more  refined  analysis  indicates  in  virtually  any 
situation  the  existence  of  parametric  forcing  terms  which  involve  both 
the  force  and  the  resulting  motion  together  and  constitute  homogeneous 
terms  with  time  dependent  coefficients.  In  fact  even  finer  analysis 
will  show  that  these  terms  in  turn  owe  their  origins  to  nonlinear 
coupling  terms  between  two  types  of  motion  where  one  of  the  types 
involved  has  been  approximated  by  a chosen  time  function.  Thus  the 
classical  problem  of  parametric  excitation  (Belaiev's  problem)  of  a 
beam  under  periodic  axial  excitation  is,  at  a more  precise  level,  a 
question  of  the  nonlinear  interaction  of  axial  and  bending  motions: 
at  a lower  level,  where  the  axial  displacement  is  taken  as  a known 
function  of  time  resulting  from  the  given  excitation,  it  becomes  a 
problem  of  bending  motion  with  parametric  excitation. 

Parametric  excitation  terms  behave  like  quadratic  nonlinearities 
in  many  ways,  this  feature  is  brought  out  in  the  report.  In  the 
heirarchy  of  'small  terms'  quadratic  nonlinearities  and  parametric 
terms  have  about  the  same  placing  although  of  course  the  parametric 
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terms  are  basically  linear.  If  the  constant  coefficient  structural 
model  with  ordinary  forcing  is  regarded  as  the  most  primitive  form 
then  the  next  stage  of  refinement  is  that  model  with  quadratic 
nonlinearities  and  parametric  forcing  included.  It  is  about  such  a 
level  of  modelling  and  the  phenomena  which  it  introduces  under  certain 
types  of  excitation  that  this  report  is  concerned. 

Excitation  itself  can  be  subdivided  into  two  main  categories, 
stochastic  and  deterministic.  These  in  turn  can  be  reduced,  stochastic 
to,  among  others,  stationary  wide  band,  stationary  narrow  band,  and 
non-stationary,  while  deterministic  excitation  can  be  monofrequency 
sinusoidal,  multifrequency  periodic  (square  wave  for  instance),  multi- 
frequency quasi-periodic  (such  as  ^ f^cosfl^t)  and  non-periodic  or 

n 

transient . 

If  attention  is  restricted  to  deterministic  excitation  then 
both  the  ordinary  and  the  parametric  excitation  can  exist  in  any  of 
the  above  categories.  There  is  clearly  a wide  range  of  possible 
forced  vibration  problems  and  very  little  is  known  about  any  of  them. 

The  present  investigation  considers  both  the  ordinary  and  parametric 
excitation  to  be  of  the  simplest  type,  monofrequency  sinusoidal,  and 
also  considers  them  to  have  the  same  frequency.  The  range  of 
phenomena  arising  if  they  have  different  frequencies  or  if  one  or  other 
or  both  of  them  is  a multifrequency  input  is  very  greatly  extended. 

In  the  analysis  which  follows,  the  linear  normal  modes  of  the 
structure  have  been  used  as  coordinates.  For  small  nonlinearities 
this  seems  a reasonable  choice,  for  highly  nonlinear  problems  some 
sort  of  nonlinear  normal  mode  might  be  more  appropriate. 

The  report  discusses  two,  three  and  four  mode  models  but  it 
should  be  understood  that  the  analysis  is  not  restricted  to  this 
number  of  degrees  of  freedom  in  the  whole  structure.  These  numbers 
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represent  the  number  of  interacting  modes  in  the  phenomena  being 
described.  Thus  the  equations  relating  to  the  "two-mode  model"  would 
apply  to  say  the  10th  and  17th  modes  of  a complex  structure  if  these 
modes  had  natural  frequencies  in  or  near  the  necessary  2:1  ratio. 

1.2  SPECIFIC 

As  mentioned  in  the  previous  section,  this  investigation  is 
concerned  with  the  resonant  behaviour  of  general  multimode  structural 
systems,  inclusive  of  quadratic  inertial  nonlinearities  and  linear 
parametric  terms,  subjected  to  monofrequency  sinusoidal  excitation. 

A formal  derivation  of  the  relevant  general  equations  of  motion  is 
given  by  Barr  5 Nelson  [1]  who  consider  purely  autoparametric 
interactions  in  structures  undergoing  harmonic  excitation  near  a natural 
frequency,  and  this  paper  forms  the  basis  for  the  current  work.  However, 
the  approach  presented  here  is  more  general  than  that  given  in  [1]  as 
parametric  terms  are  included  in  the  analysis  and  consideration  is 
given  to  other  possible  external'  resonance  phenomena  that  arise  quite 
apart  from  the  direct  forcing  of  one  mode. 

Qualitative  experimental  evidence  illustrating  the  resonance 
behaviour  exhibited  by  simple  structures  having  appropriate  internal^ 
resonance  conditions  is  presented  in  Section  2.  The  experimental 
approach  is  considered  valuable  in  that  it  demonstrates  the  physical 
existence  of  the  phenomena  and  provides  guidelines  upon  which  an 
appropriate  theory  can  be  formed.  Experimental  evidence  showing 
simultaneous  resonant  behaviour  in  up  to  four  modes  is  given. 

Some  integer  relationship  between  the  excitation  frequency  and  the 
natural  frequencies. 

Some  integer  relationship  between  the  natural  frequencies  themselves. 
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The  general  theoretical  approach  outlined  in  Section  3 is  the 
result  of  development  through  application  to  idealised  mathematical 
models,  which  are  based  on  the  experimental  structures.  The  analysis 
of  general  two  and  three  mode  structural  systems  gives  sufficient 
background  for  the  theory  to  the  expanded  to  any  multimode  system 
having  quadratic  inertial  nonlinearities  and  parametric  terms. 

Certain  modifications  to  the  basic  'Struble'  [2]  solution  technique 
have  been  found  to  be  necessary  in  order  to  maintain  reasonable 
accuracy  in  the  solution  and  to  avoid  problems  that  arise  when  both 
external  and  internal  detuning  is  considered.  The  form  of  the 
generating  solution  in  the  Struble  analysis  is  modified  by  selecting 
an  unknown  response  frequency  for  the  fundamental  term  in  the  solution, 
rather  than  the  linear  natural  frequency,  cf  [1],  although  this 
unknown  frequency  is  taken  to  be  close  to  the  natural  frequency.  In 
certain  cases  a D.C.  term  (zero  frequency)  and  a forced  term  at  the 
excitation  frequency  are  also  included  in  the  generating  solution. 

In  the  first  approximation  used  the  technique  is  essentially  a limited 
harmonic  balance  with  selection  and  retention  of  the  more  important 
nonlinear  terms. 

The  results  of  the  analysis  are  the  'variational'  equations 
which  are  nonlinear  first  order  differential  equations  describing 
the  slowly  varying  amplitudes  and  phases  of  the  solution.  For  each 
resonance  situation  in  a given  structure  there  is  a corresponding  set 
of  variational  equations,  and  all  the  possible  sets  for  two  and  three 
mode  structures  are  presented.  From  these  the  equations  for  any 
multimode  system  can  be  deduced.  Generally  the  variational  equations 
can  only  be  solved  numerically  for  a given  structural  system  and 
examples  of  the  application  of  the  theory  to  structures  having  two 
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and  four  basic  modes  are  given.  The  results  are  compared  with  the 
direct  integration  of  the  full  equations  of  motion. 

Further  analysis  of  a simplified  form  of  the  variational  equations 
in  the  'hard'  forcing  cases,  i.e.  when  no  mode  is  being  directly  excited, 
yields  regions  in  the  excitation  amplitude- frequency  plane  within 
which  the  system  will  exhibit  resonant  behaviour.  Recourse  has  to 
be  made  to  the  integration  of  the  complete  variational  equations  in 
order  to  study  the  actual  behaviour  within  the  regions  predicted. 


I 
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2. 


EXPERIMENTAL  WORK 


2.1  EXPERIMENTAL  EQUIPMENT 

Figure  2.1.1  is  a schematic  layout  of  the  experimental  equipment. 
The  model  structure  under  investigation  is  clamped  to  the  armature  of 
a medium-sized  electrodynamic  shaker  (1500  lb  maximum  vector  thrust). 

The  shaker  is  driven  by  an  accurate  sinusoidal  signal  generator 
operating  through  a power  amplifier.  Excitation  intensity  is  monitored 
by  an  accelerometer  mounted  on  the  moving  armature.  The  servo-control 
sweep  oscillator  can  be  used  in  place  of  the  signal  generator  to 
perform  programmed  frequency  sweeps  with  predetermined  excitation  levels, 
but  it  is  also  useful  in  that  its  inbuilt  signal  conditioning  circuits 
provide  a display  of  the  actual  excitation  level,  giving  the  peak  level 
of  acceleration,  velocity  or  displacement. 

The  model  structures  are  small  in  comparison  with  the  size  of 
the  shaker.  This  is  to  minimise  the  effect  of  any  feedback  from  a 
resonating  structure  causing  distortion  of  the  excitation  waveform. 

The  response  of  the  models  themselves  is  monitored  by  strategically 
placed  accelerometers  and  strain  gauges.  The  signals  from  these  sensors 
are  conditioned  through  standard  bridge  and  preamplifier  circuits  to 
give  voltage  signals  that  can  be  directly  related  to  the  behaviour  of 
the  structure.  These  resulting  signals  are  manipulated  by  simple 
analog  summing  and  scaling  devices  in  such  a way  as  to  give  signals 
representative  of  the  'modes'  of  the  system.  In  free  vibration  each 
of  these  signals  has  just  one  frequency  component,  that  is  the  natural 
frequency  of  the  mode  concerned.  In  forced  vibration  the  signals  can 
have  two  frequency  components,  one  near  the  modal  frequency  and  the 
other  at  the  excitation  frequency.  Where  required  the  active  filter 
is  used  to  remove  the  forced  response  component. 
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The  behaviour  of  the  modes  of  the  system  are  then  displayed  on 
the  oscilloscope  and/or  recorded  on  the  ultraviolet  recorder.  A.C. 
voltmeters  are  used  in  the  production  of  'steady  state'  response  curves. 
Response  frequency  measurement  is  achieved  by  the  production  of  a 
stationary  Lissajou's  figure  on  an  oscilloscope  (involving  the  use  of 
a further  accurate  signal  generator) . 

2.2  TWO  MODE  MODEL 

Figure  2.2.1  shows  a photograph  and  drawing  of  the  two  mode  model. 
The  construction  is  such  that  the  effective  lengths  of  both  steel  strip 
beams  are  adjustable.  By  careful  adjustment  of  the  masses  and  the 
lengths,  the  two  lowest  natural  frequencies  are  set  in  a 2-1  ratio. 
Figure  2.2.2  shows  the  natural  frequencies  and  the  type  of  mode  shapes 
obtained.  Strain  gauge  pairs  sited  at  the  roots  of  each  beam  provide 
signals  which  over  a wide  range  are  found  to  be  linearly  related  to  the 
end  displacement  of  the  beam  under  static  end  loading.  These  signals 
are  directly  related  to  the  generalised  coordinates  to  be  used  in  the 
theoretical  analysis  of  the  structure. 


The  separation  of  the  modes  from  these  two  signals  is  achieved 
by  their  manipulation  in  a purpose-built  analog  unit.  Figure  2.2.3 
shows  the  schematic  layout. 


SC 


Two  Mode  Experimental  Model 


Figure  2«2»2  Natural  Frequenciee  and  Mode  Shapes 


Dimensions  in  mm 
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In  effect  the  analog  devices  perform  the  same  transformation  to  normal 
coordinates  and  P2  from  the  generalised  coordinates  and  W2  as 
does  the  modal  matrix  in  the  corresponding  linear  theory. 

Figure  2.2.4  is  a copy  of  an  ultraviolet  trace  showing  the 
behaviour  of  the  model  when  the  excitation  frequency  is  close  to  the 
lower  natural  frequency.  Starting  from  rest,  growth  in  the  lower 
mode  is  rapidly  followed  by  growth  of  the  higher  mode.  Figure  2.2.5 
shows  the  steady  state  response  curves  obtained  in  the  neighbourhood 
of  this  external  resonance  condition  whilst  maintaining  the  excitation 
level  at  0.1  mm  peak  and  manually  sweeping  up  and  down  the  frequency 
range  from  8.5  to  9.1  Hz.  The  ordinary  linear  response  curve  of  mode  1 
is  replaced  by  the  twin-horned  response  curve,  accompanied  by  a similar 
curve  for  mode  2. 

Figure  2.2.6  is  of  the  ultraviolet  record  taken  when  the 
excitation  frequency  is  in  the  neighbourhood  of  the  higher  natural 
frequency.  Growth  in  the  higher  mode  leads  to  growth  in  the  lower 
mode  through  the  internal  resonance  condition.  Figure  2.2.7  shows 
the  steady  state  response  curves  obtained  in  the  neighbourhood  of  this 
resonance.  In  this  case  discontinuities  occur  in  the  stationary 
response  curves.  For  increasing  frequency  mode  2 follows  the  linear 
resonance  curve  to  point  a.  At  this  point  mode  1 jumps  from  zero  to 
a high  amplitude.  Mode  2 then  follows  the  curve  abc  and  at  c jumps 
back  to  the  declining  linear  resonance  curve.  For  decreasing 
frequency  mode  2 follows  the  curve  dbae,  having  discontinuities  at  b 
and  e where  mode  1 jumps  in  and  collapses  respectively.  It  should 
be  noted  that  only  discontinuities  in  slope  occur  in  mode  2 at  points 
a and  b.  This  'overlapping'  of  the  response  curves  for  increasing 
and  decreasing  frequencies  is  similar  to  that  illustrated  by  Haxton  [3] 
in  connection  with  the  Autoparametric  Vibration  Absorber.  The  slight 
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skewing  of  response  curves  about  the  tuned  frequency  is  attributable 
to  the  two  modes  having  natural  frequencies  that  are  not  quite  in  a 
2:1  ratio,  that  is  the  internal  resonance  condition  is  not  precisely 
met. 

Figure  2.2.8  shows  the  time  history  when  excitation  is  at  the 
combination  frequency,  + 0)2-  Excitation  level  is  high  in  order  to 
achieve  large  amplitudes  for  nonlinear  effects  to  become  dominant, 
and  the  record  starts  with  the  ordinary  forced  response,  at  the  forcing 
frequency,  of  each  mode.  There  is  then  a transformation  in  response 
as  parametric  and  nonlinear  effects  take  over.  It  should  be  noted 
that  in  this  case,  by  virtue  of  the  internal  resonance  condition,  the 
forcing  frequency  is  not  only  near  the  sum  of  the  two  natural 
frequencies,  but  is  also  three  times  the  lower  natural  frequency. 

Also  for  this  case  moderate  internal  detuning  of  the  system  was  found 
to  result  in  desynchronisation  resulting  in  a continuous  beating  between 
the  modes. 

Figure  2.2.9  shows  the  behaviour  of  the  system  when  excited 
near  twice  the  higher  natural  frequency.  Again  the  excitation  amplitude 
is  high.  The  higher  mode  is  excited  'parametrically ', leading  to  the 


excitation  of  the  lower  mode  through  the  internal  resonance  condition. 

The  result  is  a continuous  interchange  of  energy  between  the  modes. 

In  this  case  the  external  tuning  v’s  found  to  be  critical  and  the  quasi 
steady  state  shown  was  only  obtained  at  a single  tuned  frequency. 

External  detuning  over  a narrow  band  leads  to  similar  behaviour  but  of 
a more  erratic  nature.  The  'beating'  between  the  modes  becomes  irregular. 
Again,  as  a consequence  of  the  internal  resonance  condition  the 
excitation  frequency  is  not  only  twice  the  higher  mode's  frequency 
but  also  four  times  that  of  the  lower  mode. 
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2.3  FOUR  MODE  MODEL 

Figure  2.3.1  shows  the  experimental  configuration  of  the  four 
mode  model.  It  resembles  the  'T'  tail  of  an  aircraft  structure.  Again 
as  all  the  beam  lengths  are  adjustable  it  is  possible  to  set  up  different 
modal  frequencies.  In  a first  approximation  the  continous  structure 
can  be  idealised  as  a lumped  mass  system  having  'weightless'  beams 
and  four  concentrated  end  masses.  By  considering  only  ordinary  bending 
of  the  elastic  members  to  occur  then  the  system  is  described  by  four 
generalised  displacements,  and  hence  it  has  four  natural  frequencies 
with  corresponding  mode  shapes.  These  can  be  divided  up  into  modes 
involving  symmetric  and  antisymmetric  'tailplane'  motions,  as  shown 
in  Figure  2.3.2. 

Having  set  up  an  arbitrary  ' fin/tailplane'  configuration  and 
measured  the  antisymmetric  frequencies,  the  fuselage  length  is  adjusted 
to  make  the  highest  symmetric  mode  frequency  as  near  as  possible  equal 
to  the  sum  of  the  two  antisymmetric  mode  frequencies.  This  constitutes 
a combination  internal  resonance  condition  between  three  of  the  modes  of 
the  structure.  Figure  2.3.2  shows  the  frequencies  and  approximate 
mode  shapes  obtained. 
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The  behaviour  of  the  model  is  monitored  through  the  accelerometer 
at  the  base  of  the  ’fin'  and  through  strain  gauge  pairs  sited  at  the 
root  of  each  beam  forming  the  'T'  tail.  The  signals  from  these 
represent  the  generalised  coordinates  to  be  used  in  the  theoretical 
analysis.  Symmetric  and  antisymmetric  tailplane  motions  are 
monitored  by  summing  and  differencing  the  two  tailplane  signals.  It 
was  found  that  the  accelerometer  signal  and  the  summed  tailplane 
signal  were  reasonably  representative  of  the  two  symmetric  modes, 

Pj  and  P2  respectively.  However  the  fin  and  antisymmetric  tail  signals 
required  further  manipulation  to  uncouple  them  and  to  give  signals 
representative  of  the  antisymmetric  modes  P^  and  P^.  Figure  2.3.3 
shows  the  schematic  layout  of  the  equipment. 


The  variable  resistors  are  adjusted  to  that  in  free  vibration  the 
signals  P^  and  P^  contain  the  single  frequencies  Wj  and  only. 


L 
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As  mentioned  earlier,  due  to  the  nature  of  the  geometry  of  the 
structure  it  is  possible  to  excite  only  the  symmetric  modes  directly  by 
motion  of  the  shaker  head.  Apart  from  imperfections  of  construction, 
the  antisymmetric  modes  can  be  excited  by  parametric  or  autoparametric 
action  only.  Figure  2.3.4  shows  the  ultraviolet  trace  taken  when 
excitation  is  near  the  highest  natural  frequency.  The  growth  of  the 
highest  mode,  which  by  linear  theory  would  be  expected  to  predominate, 
rapidly  degenerates  as  the  first  and  third  modes  are  excited  through 
the  internal  resonance  condition.  The  second  mode  is  also  seen  to 
exhibit  small  amplitude  oscillations  at  about  half  the  excitation 
frequency  (note  also  by  coincidence  at  about  twice  that  mode's  natural 
frequency).  The  two  antisymmetric  modes  are  effectively  acting  as  an 
absorber  on  the  higher  symmetric  mode  in  a similar  fashion  to  the 
Autoparametric  Vibration  Absorber  (A.V.A)  [3].  In  this  case  better 
frequency  demultiplication  is  achieved  by  virtue  of  the  combination 
internal  resonance  condition.  [In  the  A.V.A.  a simple  two  mode  system 
having  2:1  frequency  ratio  is  considered,  similar  to  the  system 
described  in  the  previous  section].  Figures  2.3.5  (a)  and  (b)  show 
the  qualitative  steady  state  response  curves  for  modes  1,  3 and  4. 

The  former  is  the  frequency  increasing  case,  whilst  the  latter  is  the 
reverse,  the  excitation  amplitude  being  held  constant  at  0.05  mm.  As 
was  found  with  the  corresponding  two  mode  case,  jump  phenomena  occur 
and  there  is  a substantial  skewing  between  the  two  graphs,  implying 
imperfect  internal  tuning.  The  system  exhibited  unsteady  behaviour 
especially  on  the  left  hand  parts  of  the  decreasing  curves,  and  the 
points  plotted  here  are  of  the  visual  averages  taken. 

Although  the  curves  presented  here  are  qualitative,  it  is  useful 
to  have  some  idea  of  the  relative  magnitudes  of  the  responses  of  the 
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modes  involved.  By  relating  the  response  levels  of  the  modes  (mV) 
to  an  equivalent  actual  response  amplitude  of  some  part  of  the 
structure,  some  idea  of  the  relative  magnitudes  of  the  responses  can 
be  achieved.  The  antisymmetric  modes  (1  and  3)  can  be  related  to  peak 
amplitude  of  the  fin  bending  displacement  (at  the  mass),  whereas  the 
symmetric  modes  can  be  related  to  the  peak  amplitude  of  the  'fuselage' 
displacement  (at  the  accelerometer).  Hence,  as  an  approximate 
calibration: 

6 mm  total  fin  bending  excursion  corresponds  to  85  mV  in  mode  1 

and  100  mV  in  mode  3 

1.5  mm  total  fuselage  bending  excursion  corresponds  to  30  mV  in  Mode  4. 

Surprisingly  perhaps  it  was  found  that  excitation  of  high  intensity 
at  twice  the  first  or  third  modes'  natural  frequencies  did  not  produce 
any  substantial  'parametric'  resonant  behaviour  in  either  of  these  modes. 

The  only  other  resonance  behaviour  exhibited  by  the  model  apart  from 
the  ordinary  linear  forced  resonance  of  the  lower  symmetric  mode  is 
when  excitation  is  of  high  intensity  at  the  summed  frequency  of  the 
two  symmetric  modes.  This  constitutes  a combination  external  resonance 
condition.  Figure  2.3.6  is  a photograph  of  a typical  ultraviolet  trace 
taken  of  this  phenomenon.  The  second  and  fourth  modes  are  simultaneously 
excited  parametrically,  and  the  growth  in  the  fourth  mode  leads  to 
the  eventual  participation  of  the  first  and  third  modes  through  internal 
resonance.  This  is  a remarkable  situation  in  which  excitation  at  an 
apparently  harmless  frequency,  far  removed  from  any  natural  frequency, 
results  in  large  amplitude  responses  of  all  the  modes  of  the  structure. 

Figure  2.3.7  gives  a more  detailed  picture  of  the  'quasi-steady  state' 
trace  obtained  for  this  situation. 

Again  the  external  tuning  is  critical  and  slight  detuning  about 
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the  central  frequency  causes  erratic  behaviour  of  a cyclic  nature. 


Typically,  after  initial  growth  and  interaction  similar  to  the  tuned 
case  the  'beating'  behaviour  degenerates  and  decays,  whereupon  the 
process  repeats  itself. 

2.4  RESPONSE  FREQUENCIES 

For  the  purposes  of  the  development  of  the  theoretical  analysis 
it  was  required  to  determine,  where  possible,  the  actual  response 
frequencies  of  each  mode  at  each  external  resonance  condition,  and  to 
find  the  relationship  of  these  frequencies  to  the  excitation  frequency. 
The  frequencies  could  only  be  accurately  measured  for  the  cases  in 
which  steady  state  is  attained,  the  results  of  which  may  be  summarised 
as  follows. 

For  the  two  mode  system  under  excitation  near  a natural  frequency, 
the  response  of  that  mode  is  always  at  the  excitation  frequency.  Also 
in  all  cases  considered  the  higher  mode  always  responds  at  exactly 
twice  the  frequency  of  the  lower  mode,  irrespective  of  internal  or 
external  detuning.  In  the  combination  case  the  response  frequencies 
sum  exactly  to  the  excitation  frequency,  again  irrespective  of  the 
detuning  and  excitation  amplitude.  Thus  the  response  frequencies 
are  'captured'  by  the  excitation  frequency  such  that  a linear  relation- 
ship with  exactly  integer  coefficients  exists  between  them. 

In  the  four  mode  system  it  was  possible  to  collect  the  data  only 
for  the  case  of  direct  excitation  of  the  highest  mode.  Again  the 
forced  mode's  frequency  is  at  that  of  the  excitation,  and  the  first 
and  third  modes'  frequencies  sum  to  the  excitation  frequency.  However 
as  there  are  now  three  modes  involved  in  the  internal  relationship, 
the  response  frequencies  of  the  first  and  third  modes  are  not  fixed 
upon  specification  of  an  excitatioh  frequency,  the  only  condition 
being  that  they  sum  to  it.  Experimental  evidence  suggests  that  a 
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first  approximation  for  the  frequencies  can  be  obtained  by  assuming 
that  any  detuning  is  shared  between  the  modes  in  proportion  to  their 
natural  frequencies.  Thus  in  this  case,  the  appropriate  value 


for  (the  response  frequency  of  the  first  mode)  would  be  given  by 
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3.  THEORY 

3.1  EQUATIONS  OF  MOTION 

The  general  form  of  the  equation  of  motion  of  the  rth  mode  of 
a discrete  structure,  retaining  quadratic  inertial  nonlinearities 
and  linear  parametric  terms  is 

Pr  ^ - «[Ki.jPjCOsQ^t  . ^ "'rijPiPj  " ‘"ArPr^ 

3.1.1 

where  summation  on  the  repeated  indices  i,  j only  is  implied.  Examples 
of  the  derivation  of  these  equations  for  specific  structural  configurations 
are  given  in  the  next  section.  is  the  direct  forcing  term,  the 

parametric  term,  both  at  frequencies  while  j » 

quadratic  inertial  nonlinear  terms,  e.  being  a small  parameter.  As  in 
[1]  it  is  assumed  for  simplicity  that  the  damping  is  viscous  and  that 
there  is  no  damping  coupling  between  the  modes. 

3.2  SOLUTION  TECHNIQUE 

A solution  is  taken  from  the  form 

Pr  = Aj.o(t)  + A^(t)cos(fl^t  + (J)^(t))  + f^cos(fljt  + + 6a^(t)  + £^b^(t)  .. 

3.2.1 

The  first  three  terms  are  the  principle  part  of  the  solution,  while 

the  additive  terms  in  powers  of  e allow  for  a perturbational  treatment. 

The  first  term  is  a D.C.  or  rectified  component,  the  second  a component 

with  slowly  varying  amplitude  and  phase  at  frequency  , the 

t h 

response  frequency  of  the  r mode,  as  yet  unknown.  The  third  term  is 
the  component  of  forced  response  of  the  mode  at  the  forcing  frequency. 

The  linear  undamped  response 
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factor  on  and  the  phase  ip^,  can  be  incorporated  but  this  is 
not  usually  necessary  as  the  component  is  brought  in  only  for 
excitation  frequencies  reasonably  well  removed  from  ordinary 
resonance  and  at  such  points  the  damping  plays  little  part  in  the 
response. 

The  solution  form  3.2.1  requires  certain  modifications  when  the 
excitation  frequency  is  at  or  near  a natural  frequency  of  the  system. 
When  exciting  near  the  natural  frequency,  = Qj  (the  mode  responds 
at  the  forcing  frequency)  and  there  is  duplication  in  the  response 
components.  Hence  the  linear  forced  response  of  the  mode  is  not 
incorporated  in  the  solution.  In  such  a situation,  the  excitation 
will  generally  be  'soft'.  That  is,  very  small  excitation  amplitudes 
are  required  to  generate  large  responses.  Thus  the  Fj.'s  can  be  taken 
to  be  of  order  e,  or  less,  and  so  in  a first  approximation  it  is 
reasonable  to  neglect  the  linear  forced  response  of  all  the  other  modes 
as  well.  Also  it  has  been  found  in  such  cases  that  the  D.C.  terms 
have  little  influence  on  the  accuracy  of  solution.  Hence  the  solution 
from  when  the  excitation  is  'soft'  reduces  to: 

Pr  = A^(t)cos(n^t  + (})^(t))  + ca^(t)  + 6^b^(t)  + ...  3.2.3 

which  is  similar  to  that  used  by  Barr  et  al  [1],  [3],  [4],  [5]. 

However,  when  the  excitation  is  'hard',  i.e.  when  it  is  away 
from  a natural  frequency  and  at  an  amplitude  that  can  no  longer  be 
taken  to  be  small,  it  has  been  found  that  both  the  linear  forced 
response  and  the  D.C.  components  are  essential  to  the  analysis  if 
reasonable  accuracy  of  this  first  order  solution  is  to  be  maintained. 

Taking  the  most  general  case  ('hard'  forcing),  the  solution 
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procedure  begins  by  the  substitution  of  3.2.1  into  3.1.1.  Restriction 
of  the  analysis  to  one  involving  terms  of  the  first  order  in  s results  in: 


Aro  [^r  - ^ ]cos (5^^t  + ())^)  - [A^$^  + 2A^(S2^  + ^^)]sin(fi^t  + 0^)  - 

- f^Q|cos(i2jt  + + ea^  + + A^cos  + 4,^)  + f^cos(f2^t  + ip^)  + 

+ ea^J  = - cjp  [A^^  + A^cos(J^^t  + (j)^)  + f^cos(fi^t  + + ea^]  + 

+ F cosf2,.t  - e{K  [A.  + A.cosfri.t  + A . ) + f.cos(f2^t  + iii.)]cosfi^t  + 
r t ^ rj  ^ JO  j '■  j j '■  f ^j-'-'  f 

+ I . . [A.  + A.  cos(S2.  t + 4 . ) + f . cos  (f2^t  + il/ . ) ] • [A  • + (A.  - A . (fi . + 4>  ■ ) ^)  • 

rij 10  1^1  1 '■  f ^1'^-'  JO  '■  J J J ^j-* 

.cos(fi.t  + (p.)  - (A.$.  + 2A.(fi.  + d) . ) }sin  .t  + c]).)  - f.fi|cos(Q^t  + ij;.)]  + 

+ m^ij[A.^  + A^cos(Q^t  + 4)^  - + ^^)sin(n^t  + 4).)  - + 4;^)]. 

.[A.  + A.cos(fi.t  + 4>-)  - A.ffi.  + 4).)sin(J2.t  + 0.)  - f.f2„sin(fi^t  + 4;.)]  + 

i JO  J V J Tji  ji  J Yji  1 J Yji  J f >.  f ryj 

+ 0)  d [A  + A cos(fi  t + 4>)-A(f2  + 4)  )sin(fi  t + 4)^)  - f ^^sinCfi^t  s+  4^^)  ] } 

r rr  ^ ro  r r r r r r r r r r ^ 


3.2.4 


In  order  to  make  the  analysis  manageable  at  this  stage  these 
equations  are  not  treated  precisely,  and  the  slowly  varying  accelerations 
and  products  of  velocities  are  neglected  (A^,  0^,  A^^^  etc.). 

Similarly,  the  slowly  varying  velocities  of  order  e are  also  neglected. 

The  analysis  continues  by  considering  terms  of  equal  powers  in  e. 
The  zero  order  terms  form  the  so  called  'variational  equations'.  Terms 
which  cause  resonance  in  the  first  order  perturbational  equation  (c*) 
are  removed  and  considered  along  with  the  variational  equations.  Hence 
the  resulting  'variational'  equations  contain  the  essential  elements 
that  control  the  behaviour  of  the  system.  The  D.C.  terms  and  the  cos 
and  sin  terms  in  (fi^t  + 4*^)  are  considered  separately.  In  the  first 
approximation  the  technique  is  essentially  a limited  harmonic  balance 
with  preferential  selection  and  retention  of  the  more  important 


nonlinear  terms. 


From  3.2.4  the  fundamental  variational  equations  are: 


A wi  = 0 


(D.C.  terms) 


3.2.5 


- wpA^  (cosCn^t  + (t)^)  terms)  (3.2.6 


■2A  Q =0 
r r 


(sin(l^  t + (})  ) terms 
r r 


3.2.7 


and  the  first  order  perturbational  equations  are: 


a + u)^a  = - { 
r r r 


K .[A.  + A.cos(fl.t  + ({).)  + f.cos((l„t  + ij; . ) Icosflj^t  + 

rj  ^ JO  j '■  j J ^ f ■“  f 


+ + A^cos(fi^t  + i))^  + f^cos(fi^t  + i|;^)  ] . [-AjOjCOS  (Hjt  + (|)j)  - f^ 


l^cos(fi^t  + ip^)]  + [-A^fi^sin(J2^t  + 4>^)  - f^^^^sin (J)^t  + t);.)] 


[-A^Jl.sin((ljt  + 4).)  - f^fl^sin(Q^t  + 4/^)]  + U3^d^^ [-A^fi^sin (fi^t 


+ 4>j.)  - f^n^sin(Jl^t  + 4;^]} 


3.2.8 


It  must  be  noted  that  when  the  forcing  is  'soft',  the  term 
(order  €)  should  also  be  included  in  this  equation  as  it  will  be  a 
resonant  term  for  the  mode  that  is  being  directly  excited.  However  in 
that  case  the  D.C.  terms  (A^^) , and  the  linear  forced  response  terms  (f^) 
will  be  omitted. 

The  first  damping  term  in  3.2.8  (second  last  term  on  RHS  of  the 
equation)  is  clearly  resonant  and  must  be  removed  to  the  corresponding 
variational  equation.  D.C.  terms  arise  where  there  is  a product  of  two 
trigonometric  functions  having  the  same  frequency.  For  example 
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cos(fi^t  + iJ^j)cosn^t  expands  to  J[cos(2fij.t  + + cos^Jj  ] , a harmonic 

part  at  2fi^  and  a D.C,  part  (ij/j  constant  for  a given  excitation 
frequency).  All  such  D.C.  terms  are  removed  to  the  D.C.  variational 
equations.  Other  terms  can  become  resonant  if  there  are  certain 
relationships  between  the  forcing  frequency  and  the  responding 
frequencies,  and  amongst  the  responding  (near  natural)  frequencies 
themselves . 


3.3  EX-VERNAL  AND  INTERNAL  RESONANCE 

On  examination,  equations  3.2.8  yield  the  frequency  relationships 

that  may  cause  resonance  in  the  system.  Consider  the  parametric  term, 

(•••)>  then  as  cosfl^tcos (iljt  + (jjj)  expands  to  |[cos(il^t  + fl^t  + ({>j)  + 

cos(fl^t  - iljt  - c|)j)],  this  term  can  be  resonant  in  the  j ' mode  if 

= 2fij  or  in  the  i^^^  mode  if  ± . These  frequency  relations 

are  between  the  excitation  frequency  and  either  one  or  two  of  the 

modal  response  frequencies  and  hence  are  termed  'external'  r .sonance 

conditions.  These  are  in  fact  the  ordinary  parametric  and  combination 

parametric  resonance  conditions.  The  terms  involving  products  of  the 

excitation  and  response  frequencies  in  the  H . . and  m . . 'nonlinear' 

^ rij  rij 

brackets  of  3.3.8  also  yield  the  same  conditions.  Given  the  external 

resonance  condition  these  are  the  types  of  terms  that  may  cause 

external  excitation  of  the  modes  involved  in  it. 

The  nonlinear  terms,  however,  predict  a further  type  of  resonance 

condition.  Terms  involving  products  such  as  cos(il^t  + ())^)cos(njt  + (Jj^) 

produce  fundamental  harmonics  of  ± Jl^).  These  can  be  resonant  in 

the  mode  if  Q.  = 2^-  or  in  the  r'"^  mode  if  J)  = fl.  ± Q-.  These 
1 j r 1 j 

conditions  are  termed  'internal'  resonance  conditions  as  they  are 
independent  of  the  external  excitation  but  depend  upon  the  natural 
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frequencies  of  the  system.  Given  an  internal  resonance  condition, 
the  nonlinear  terms  may  be  such  as  to  allow  a mode  that  is  already 
excited  to  excite  another  of  the  modes  involved  in  the  'internal' 
relationship  (otherwise  known  as  'autoparametric ' interaction). 


I 

To  sum  up,  the  possible  resonance  conditions  are:-  i 


(direct  'soft'  excitation) 
('parametric'  excitation) 
('combination'  parametric  excitation) 


external 

resonance 


n. 

1 

Q 


r 


2(1. 

J 

(1.  + (1 . 
1 J 


(implying  = 2Uj) 
(implying  co^  = + to^) 


internal 

resonance 


Thus,  to  the  order  e',  one  external  forcing  frequency  can  excite 

up  to  a maximum  of  two  modes.  It  is  possible  that  each  of  these  excited 

modes  in  turn  can,  if  the  internal  resonance  conditions  exist,  excite 

up  to  two  more  modes  ' autoparametrical ly' , and  so  on.  Experimental 

evidence  of  this  cascading  effect  has  been  presented  in  Section  2. 

Cubic  and  higher  order  nonlinearities  give  rise  to  other  external 

[6]  and  internal  [1]  resonance  conditions.  For  example  cubics  can 

give  external  conditions  such  as  (2^  = d-  + 2(1-,  (2.  + (1.  + Q,  , and 

f 1 j’  1 J k’ 

internal  conditions  such  as  (1  = 2(1  ■ + (2 • , (1-  + (1  - + (1.  . The  higher 

T 1 J 1 J K 

the  order  of  nonlinearity  the  more  modes  can  become  involved.  Here 
interest  is  confined  to  the  quadratic  nonlinearities  and  the  parametric 
terms  only. 


3.4  THE  VARIATIONAL  EQUATIONS 

At  this  stage  in  the  analysis  it  becomes  necessary  to  study 
particular  structural  systems.  Interest  is  focussed  on  systems  having 

I 
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internal  resonance  of  which  for  this  class  there  are  two  types,  as 
shown  in  the  previous  section.  They  involve  either  two  or  three  modes 
such  that 

2^2J  = or  < ^^2  ^ ^3) 

It  is  of  course  possible  for  both  these  conditions  to  hold  simultaneously, 
but  this  situation  is  less  likely  to  be  encountered  with  in  practice. 

For  a system  having  internal  resonance  between  two  modes  in  the 
formal,  = 211^,  the  external  resonance  possibilities  are: 

^2^  2H2 > ^2 

where  the  first  two  are  'soft'  forcing  cases  (note  that  211^  = 

and  = lip  . 

For  a system  with  three  modes  in  internal  resonance, 
the  external  resonance  possibilities  are: 

Sl£  = 11^,  9^,  Hj,  212^,  2112,  29^,  9.^  - 9^  + 9^,  9^  + 11^, 

where  the  first  three  are  'soft'  forcing  cases  (note  that  122  * = ^3’ 

12^  - 122  = 12j , 9^  - 12j^  = 122) 

Each  of  these  external  resonance  conditions  produces  a corresponding 
set  of  variational  equations.  In  each  case  these  are  derived  by 
removal  of  all  the  resonant  terms  from  3.3.8  into  the  appropriate 
fundamental  variational  equation,  3.3.5,  6,  7.  The  two  mode  systems 
are  examined  first. 

I 
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(1 ) Two  mode  systems 

Taking  the  first  case,  with  ^2  = 215^^,  and  neglecting 

in  the  solution  form  the  components  of  linear  forced  response  and  the 
D.C.  terms,  ('soft'  forcing,  see  equation  3.2.3),  the  resonant  terms 
for  mode  1 are  those  involving 

'^12’  ^112’  ^121’  ”’112'  "'121*  “1^11  ^1’ 

and  for  mode  2:- 


^21’  ^211’  "’211’  “2‘^22‘ 

The  resulting  variational  equations  are:- 

= (flj  - ojpA^  + F^cosc))^  - j [P^A2COs((t)2  - 0^  + N^A^A2COS (4)2  - 24)^)] 

-2A^fi^  = - y [-P^A2sin(4)2  - 4)^)  - N^Aj^A2sin(4)2  - 24>^)  - 2WjJljdj^Aj  + F^bin4i 

= (^^2  - ■ ‘*’1^  ■ N2A^cos(4)2  - 24)^)] 

-2^2^12  = - I [P2A^sin(4)2  - '('j)  - N2AjSin(4)2  - 24)^)  - 2u>^Q2<^22^2^ 

where  P^  = <^2 

^ ^"'112  ■"  '"l2p^l^2  ■ ^112^2  ■ ^121^1 
^2  = ^21 

^2  " *'^211  "'211^^1 

Similarly  for  = SI2  ^2  ” ^^1'  resulting  variational 

equations  are:- 


( 


- j [P  jAjC0s2(()j  + N jA^A^cos  (4)2  - 2(t)j)] 
-2A^S2j  = - I [PjAjSin24)j  - N jA^A^sin  (4)2  - 2cp^)  - 2^^jU)jdjjA^J 
= {^2  - F2‘^°®‘^2  ' \ [-N2A^cos(4>2  “ 24)jJ] 

-lA^^  = F^sin4)2  - y [-N2AjSin (4)2  - 24)^)  - 2Q.2^2'^22^2^ 


where  P^  = and  and  N2  are  as  in  3.4.1  3.4.2 

For  the  hard  foi'cing  cases,  + il2  2^2^,  the  equations 

are  considerably  more  complex,  as  the  linear  forced  response  and  the 
D.C.  terms  are  included  in  the  solution  form.  For  the  two  mode  system 
the  D.C.  'variational  equations'  are:- 


3.4.3 
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The  derivation  of  these  equations  is  not  affected  by  the  internal  and 
external  frequency  relationships,  and  so  they  are  the  same  for  both 
the  hard  forcing  cases. 

For  the  combination  external  resonance  case,  with 

^2  = the  slowly  varying  amplitude  and  phase  variational  equations, 

^r 

with  f = and  =0,  are: 

- wpAj  - [PjA2Cos(<()j  + ((i^)  + N^A^A^cos (4)2  - 24)^  - 


■2(^11^10  ^ ^21^20)^S^ 


-2A^fi^  = - y [P^A2sin(4)^  + - N^A^A2Sin(4i2  - 241^)  - 2a)^fi^d^^A^] 

-2A2J^2^2  " *-^2  ■ ‘^2^^2  ' I ^*^1  *^2^  " N2A^cos(4)2  - 24)^)  - 


■^‘•‘^212'^10  ^222'^20^^2^2^ 


-2^2^2  ^ ~ 2 tP2^1^^''*'‘^l  *^2^  ■ " 2‘^2^2‘^22'^2^ 

where  P^  = k^2  ' ^ii2^l‘^2  ’ ^121^1^f  ‘ ^122^2’^^2  * ""  ^'^22' 

•^2^2^f  ■"  ^"’112  '"l21^^1^2^f 

^211^1^^!  * ^f^  * ^"’211- 

fl^A  " ^">221  " "'212^  Wf 


^2  ^ '^21  ■ ^212^2^f  ■ *'221^2^1 


3.4.4 


i 

I 
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It  is  obvious  that  certain  similarities  exist  between  equations 
3.4.1,  2 and  4.  The  terms  involving  and  are  those  which  arise 
from  the  external  resonance  condition.  Terms  involving  are  those 
which  arise  from  the  internal  resonance  condition,  and  these  are 
identical  in  each  set.  In  the  hard  forcing  situation  the  parametric 
(Pj.)  type  terms  not  only  contain  the  ordinary  parametric  terms 
but  also  nonlinear  terms  arising  from  the  inclusion  of  the  linear 
forced  response  components  in  the  solution. 

The  D.C.  'variational'  equations  3.4.3  can,  by  inspection  of 
3.2.8,  be  written  in  a general  form  applicable  to  an  n mode  structure: - 


A = 
ro  r 


T -f- 
2 ^ ri  1 


(m^. , - Jf,f.Sl‘  + 


rij 


rij 


- 

1 J f 


rjj  rjj"  j 


i, j = 1,  . . . , n (no  sum  on  r) 


3.4.5 


Also  by  inspection  of  3.4.5,  the  additions  that  these  D.C.  terms  make 
to  the  slowly  varying  variational  equations  can  be  written  in  general 
form.  The  addition  to  the  right-hand-side  of  each  -2A^Q^({)^  equation 
takes  the  form 


V 

r 


n 

y el  . A.  A 

rir  10  r r 

1 = 1 


3.4.6 


Equations  3.4.5  and  6 cover  the  D.C.  effects  in  all  the  cases  to  be 
considered. 

The  variational  equations  for  the  final  two  mode  system,  = 29,.^ 
with  ^^2  ” are: 


- (jJpA^  - I [N^A^A2COs((()2  - 24)^)]  + 


-ZA^^j  = - j [-N^A^A2Sin(<t)2  - 241^  - 2f2j(D^dj ^A J 

-2A,Q2^'2  ^ ^^2  ■ ‘^2^^2  " 1 t’’2^2‘^°^^‘*’2  ' N2AjCOs(4)2  - 24>p]  + ^2 

-2A2i^2  ~ ~ \ [P2^2^^"^'^2  ' ~ ^^2‘^2‘^22^2^ 

where  'P 2 = *^22  ' ^l^*'221^f  ^212^2  ‘ *•"’221  "'212^S^f^ 

” ^2  ^**^222  * ^2^  ~ 2ni222^2^f^ 

and  N^  and  N2  are  as  in  3.4.1.  3.4.7 

In  all  of  these  systems  of  equations  (3.4.1,  2,  4 and  7),  the 
terms  arising  from  the  internal  resonance  condition,  the  terms, 
are  the  same.  If  the  system  possesses  no  such  internal  resonance 
condition  then  these  terms  simply  disappear,  leaving  equations  similar 
in  form  to  those  studies  by  Yamamoto  et  al  [6,  7]. 


(2)  Three  Mode  Systems 

There  are  nine  separate  external  resonance  cases  for  the  three 
mode  system.  Only  the  ones  to  be  used  subsequently  are  listed  here. 
All  the  other  cases  are  listed  in  Appendix  I. 


M 


1 


ih 

; 

f'v 


V 4 

K 


= S2^  ('soft'  forcing) 


■2A^S2^(Jj  =(.Qj-wj)Aj  - ^ [PjA^cosCct^  + (J)^)  + N ^A^A^cos  ((fi^  - $2  ‘ 


" ■ 2 [PiA^sinlf})^  + (J)^)  - N^A^A^sinCc))^  - <p2  ' 'J’l)  1^1 1 

-2A2n2^2  " ^^2  ■ ‘^2^^2  ' 1 [P2^1 ■"  '^2^  * 1>  ~ ^2  ~ 

~2^2^2  ~ ~ T ^*^1  * ^2^  ” ^*^3  ~ *^2  ” ” ^^2^2*^22^2^ 

■‘'^3^3^3  ^ ^^3  ■ * FjCos^j  - j [N^A^A^cos  ((l)^  - (J>2  - <}>2^)] 

■^^3^3  " ''3^^'^'^3  ■ 2 [N3A^A2Sin(02  - $2  ‘ ‘J’p  " ^^3“3‘^33^3^ 


where  P^  = 

^2  - ^21 


^1  " *^"'123  ■"  '"l32^S^3  ■ ^123^3  ■ ^132^2 
^2  " ‘^"*213  " "’231^^1^3  ■ *^231^1  ’ ^213^3 
*^3  ^ ■*^"'312  ■"  "’32p^l^2  ■ ^312^2  ‘ ^321^1 


3.4.8 


From  these  equations,  including  those  in  Appendix  I,  it  is 
possible  to  deduce  the  equations  for  any  given  multimode  system.  For 
example,  the  equations  for  a system  having  four  modes  and  resonant 
frequency  relationships  of: 


W 


with  12,  + 12_  = 12.  are; 
f 2 4 13  4 


- cjpA^  - j lN^AjA^cos(<^^  ' '>’3  ’ * ^1 


■ZAifij  = - j [-N^AjA^sinCft)^  - (t)^  - <t)j)  - ^Aj ] 


2A2fi24'2  = (^2  ■ *^^^2  ■ I [P2\‘=°^^‘5’2  * ^4^'^  * ^2 


’ 2 * *^4^  ■ ^^2“2‘^22'^2^ 


2A3V3  = ^^3  - “^3^^  - I IN3AiA4COs(<|)4  - 0)3  - 4>i)]  - V3 


= - 2 [-N3A^A^sin(<t>4  ' ‘<’3  ‘ 


2A4f24<t>4  = (^4  - (^J)A^  - j [P^A^cosCcI)^  + + N^A^A^cos ((t)^  - (f,^  - <))^)]  + 


-2^4  = ■ 2 [P4A2Cos(<J)2  + + N^A^A^sinCtt)^  - <p^  - 4>j)  - 2S^^a)4d^^A^] 


where  P.  = k..  - - f,  («._,. + 

2 24  244  4'-  4 f^  1^  214  4 241  f-^  2^  214  4 


■"  ®'242^f^  ■ ^3^*^234'^4  ‘ ^243^f^  ^'"244^4^4^f  ^"’214  "'241^^1^4^f 


(">242  " "’224)  Wf  " ^"'234  " '"243^^3^4^f 


and  P^  - k^2  - ^422^2^^2  * ‘ ^1^*'412S  ' %21^f^  ‘ ^3^‘^432^2  " 


■"  ^423^f^  ■ ^4^^442^2  ^424^P  ^"'422^2S^f  ^'”412  ”’421^^1^2^f 


^"’432  * '"423^^3^2^f  ■*■  *^'"442  "'424^^4^2^f 


and 


1 


and 


= ^"'134  * 

’"l43^(^3(^4 

0 q2  _ 

134  ^4 

f,  q} 

143  3 

•\3 

= ("’314  * 

’"341^^^l(^4  - 

^314(^5  - 

34l“l 

"-(’"413  " 

’"4  3p(^^^3 

*^413^3  ■ 

0 o2 

431  r 

3.4.9 

3.5  SOLUTION  OF  THE  VARIATIONAL  EQUATIONS 

The  standard  procedure  at  this  stage  is  to  seek  steady  state, 

• • 

or  stationary  solutions  by  setting  A and  4>  to  zero  [B3,  H2,  II,  Al,  Y2]. 


However,  experimental  evidence  shows  that  only  in  certain  cases  will 
stationary  conditions  be  attained.  More  generally  there  is  some  form 
of  continuous  'heating'  between  the  modes  involved. 

Yamamoto  and  Saito  [7],  when  considering  purely  combination 
parametric  systems  which  exhibit  exponential  growth  in  an  instability 
region,  assume  'normal'  solutions  by  taking  4>  = 0.  This  results  in  a 
relationship  giving  the  unknown  response  frequencies,  in  terms 
of  the  excitation  and  natural  frequencies,  and  In  other  words 

it  specifies  the  'detuning'  of  the  system  defining  by  how  much  the 
response  frequencies  differ  from  the  linear  natural  frequencies.  In 
the  event,  even  for  the  simple  systems  considered,  the  relationship 
is  not  trivial  and  an  approximation  is  employed  to  simplify  it  before 
the  continuation  of  analysis. 

In  an  allied  paper  Yamamoto  and  Hayashi  [6]  consider  two  mode 
systems  having  quadratic  and  cubic  elastic  nonlinearities  under 
harmonic  excitation.  The  external  resonance  cases  studied  are 

= 2f2j  + ^2  + fl2*  Each  resulting  variational  equation 

contains  only  one  trigonometric  function  (because  of  the  lack  of 

internal  resonance  and  parametric  terms),  and  so  analytical  solutions 

• • 

are  obtained  upon  setting  A and  ^ to  zero.  Again  the  detuning 
relationships  are  quite  complex,  but  for  a first  approximation 


1 


40 


The  ratio  of  response  frequency  to  excitation  frequency  is  equal  to  the 
ratio  of  natural  frequency  to  the  corresponding  'tuned'  exciting 
frequency.  In  other  words  the  'external'  detuning  is  shared  proportion- 
ately between  the  modes.  Experimental  evidence  suggests  that  this  is 
a reasonable  approximation. 

Where  there  are  both  external  and  internal  resonance  conditions 
the  situation  becomes  more  complex.  For  two  mode  systems  the  analysis 
procedure  has  already  set  the  response  frequencies  exactly.  For 
example  3.4.2  has  been  derived  by  taking  = ^2  with  = 20^.  For 
a given  excitation  frequency  the  response  frequencies  are  fixed.  In 
the  combination  case,  + ^2  but  ^2  = and  so  again  and  f^2 

are  fixed.  These  frequency  relationships  are  confirmed  experimentally. 

For  three  mode  systems  with  soft  forcing,  the  response 
frequency  of  the  directly  resonant  mode  is  fixed  at  the  excitation 
frequency,  but  the  only  other  condition  on  the  other  modes  involved 
through  the  internal  resonance  relationship  is  that  11^  + 0^= 
and  so  the  frequencies  12^  and  are  not  fixed  exactly.  In  the  hard 
forcing  combination  situation  none  of  the  frequencies  are  set 
exactly  through  the  solution  procedure  (e.g.  = 122 

+ 122  = 12j) . 
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It  has  not  been  found  possible  to  obtain  general  analytical 
solutions  to  any  of  the  sets  of  complete  variational  equations  listed 
in  Section  2.5  and  Appendix  I. 

As  they  do  not  lend  themselves  to  analytical  solution,  recourse 
has  to  be  made  to  numerical  solution  for  particular  structural 
systems.  Where  necessary,  equations  of  the  type  3.5.1  are  used  to 
give  values  for  unfixed  frequency.  Where  steady  state  conditions  are 
found  to  exist,  by  numerical  integration,  then  it  is  possible  to  compute 

frequency  response  curves  by  numerical  solution  of  the  nonlinear 

• • 

algebraic  equations  that  result  from  setting  A^  and  (J>^  to  zero.  The 
accuracy  of  the  solution  procedure  can  be  checked,  in  the  simpler 
cases  at  least,  by  direct  integration  of  the  full  equations  of  motion 
under  the  same  excitation  conditions. 

Analytically,  more  progress  can  be  achieved  if  certain 
assumptions  are  made.  In  the  soft  forcing  cases  analytical  solutions 
can  be  found  if  the  'parametric'  terms  are  neglected.  This  gives 
results  similar  to  those  presented  by  Barr  and  Nelson  [1],  these  will 
not  be  discussed  further  here.  In  the  hard  forcing  situations  the 
initial  resonant  behaviour  is  governed  by  the  'parametric'  type  terms 
in  the  variational  equations.  The  terms  arising  from  the  internal 
resonance  conditions,  or  'autoparametric ' terms  only  become  dominant 
once  growth  in  a relevant  mode  has  begun.  Hence  by  neglect  of  these 

I 

terms  it  is  possible  to  derive  expressions  which  will  give  the 
approximate  stability  boundaries  for  the  system.  Recourse  has  to  be 
made  to  the  full  variational  equations  in  order  to  study  the  resonant 


behaviour  of  each  mode  involved. 
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3.6  STABILITY  BOUNDARIES  - hard  forcing  cases 

For  vanishingly  small  A.  allowing  the  neglect  of  the  auto- 

t > J 

parametric  coupling  terms  in  the  complete  variational  equations,  there 

are  regions  in  the  excitation  amplitude/ frequency  parameter  space  for 

these  cases  inside  which  the  A.  . are  subject  to  exponential  growth 

t > J 

with  time.  In  this  section  the  boundary  of  this  region,  which  is 
effectively  the  boundary  of  the  region  inside  which  growth  and 
interaction  occur,  is  found.  The  external  resonance  conditions  under 
consideration  are: 


+ Q.  or  = 2Q. . (j  > i,  i / j) 

f J—  1 f 1 


Taking  the  first  of  these  and  neglecting  the  autoparametric  terms, 
the  variational  equations  for  the  i^^  and  j^*^  modes  of  an  n mode 
structure  are:- 


.2A.I1.$.  = - w^)A.  - j P.A.cos(4).  1 ()).)  + ^^kAo^^i 


-2A^f2^  = * 2 P^A^sin((t)j  ^ (}).)  + efi^aj.d.^A^ 

n 

-2A.fi.^.  = (12]  - a3])A.  - j P.A.cos((^.  1 0.)  * .A^^A.Il] 

-2Ajfi^  = - j P^A^sin(0^  ^ ())^)  + eJl^w^d^^A^ 

where  P.  = k.  . + ^f^I2|  - «-.^.f^I2]  + (m.  + "'ikj ^ ^ 

’j  - "ji  ^ ji  l-"jik^“f  - '^jkifk^i  ^ ^ik  ^ "jk!^  WfJ 


(1) 


(2) 


(3) 


(4) 


and  P.  = 


3.6.1 
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These  can  be  deduced  from  equations  16,  7,  and  8,  the  upper  and  lower 
signs  being  adopted  for  summed  and  difference  combination  conditions 
respectively.  The  D.C.  terms  are  given  by  3.4.5,  which  here  take  the 
reduced  from 


A,  < = 
ko  k 


^ [ic,  f + (m,  - )f  f Qi] 

2 '■  kp  p kpq  kpq  P 1 f 


where  p,q  = l,n 


3.6.2 


as  the  terms  involving  the  response  amplitudes.  A?,  are  omitted  zero 
at  the  onset  of  instability). 

Let 


A. 

1 


1 


0)? 

1 


I ^^iki\o 

k=l 


3.6.3 


Then  taking  a solution  of  the  form 


= 0 


3.6.4 


equations  3.6. 1(1)  and  (3)  yield 

3.6.5 

3.6.6 

Hence  the  amplitude  ratios  are  fixed,  implying  common  growth  rates. 

Let 


A^ 


A' 


A.P. 

J 1 
A.P. 
1 J 


and  A. A. 
1 J 


-r-  P.P.cos^  ($.  + ({>.) 
4 1 j 


A.  . 
i.J 


3.6.7 
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r ■ 


<*■ 


then  3.6. 1(2)  and  (4)  yield: 


e ^ P . P . 

A ] 2 

:r^  sin  (d) . + (b.)  = 

4 I _ ■*'1^ 


+ (2ari!.  + efl.cj.  d . . ) (2otU . + eSl.w.d..) 
1 1 1 11-' ' j J J JJ 


3.6.8 


For  P^,  Pj  non-zero,  3.6.6  and  8 give 


g A. A.  ^2  P-P. 

+^((i).d..  +a).d..)cx+  - d.  . w.  d . . - .i  I -)  = 0 

2^111  111-^  — 4i2.fi.  4^111111  4fi.fi. 

J JJ  J J JJ  J J 


having  roots 


a = 


-p  (w.  d.  . + 0)  .d . . J 

4 ''  1 11  J jj'  W 


2 A. A. 

^(w.d..  -(D.d..)^+  Vo  n 
16^111  jjj  4fi^fij 


e^P.P.nii 
1 J 


16fi.fi. 
1 J 


3.6.9 


The  circled  ± sign  does  not  refer  to  summed  and  difference  combination 
conditions.  The  critical  root  as  regards  stability  is  the  more  positive 
one  [Yl],  viz. 


e I c ^ A . A . 

a = - T ((jD.d..  + UJ.d..)  + I—  (oj.d..  - oj.d..)^  + a- 
4 '■  1 11  J jj'  16  ' 1 11  J JJ-”  4fi.fij 


e"P.P. 


16fi . fi . 

1 3 


3.6.10 


For  a > 0,  growth  will  occur  leading  to  resonant  behaviour.  Hence  the 
stability  boundary  is  defined  by 


+ e^P.P.  = 4[e2(jj.d.  .01. d.  .fi.fi.  + A. A.] 
1 J 1 11  J J1  1 J - 1 J 


3.6.11 


the  equality  being  replaced  by  > for  growth. 
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Again,  upper  and  lower  signs  are  adopted  for  summed  and  difference 
conditions  respectively.  For  a given  structure,  3.6.11  is  a function 
of  excitation  frequency  and  amplitude.  Hence  for  a given  resonance 
condition  it  is  possible  to  compute  stability  diagrams  in  the  excitation 
amplitude- frequency  plane. 

The  governing  equations  for  the  simpler  case,  when  are: 


-ZA.Jl.cf).  = (fit 
111  1 


2 2,,  — P . A.  cos2(f) . + y e£..-A,  A.fi? 

f-coT)A.-2  11  ^1  ikikoii 


-2A.Q.  = - ^ P.A.sin2(})-  + efi.w.d..A. 

11  2ii  1 11111 


where  P.  = k..  . ' ^iik^k^f  " ^^ik  " ^ki^fkS“fJ 


3.6.12 


Repetition  of  the  preceding  analysis  yields  the  growth  condition: 


e^P.  ^ 4 [e  ^to?d?  .SI?  + A?] 
1 1 11  1 1^ 


3.6.13 


It  must  be  noted  that  the  D.C. terms  are  functions  of  the  response 
amplitudes,  A^,  and  when  approaching  the  stability  boundaries  from 
within  the  unstable  region  the  stability  boundaries  will  be  affected 
accordingly.  The  anlysis  presented  here  is  limited  to  approaching 
the  boundaries  from  the  'outside'  only. 


w 
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4.  STRUCTURAL  EXAMPLES 

4.1  INTRODUCTION 

This  section  presents  a brief  report  of  the  analysis  and  selected 
numerical  results  from  the  theoretical  models  which  are  based  on  the 
experimental  models  of  Section  2.  At  each  resonance  situation  the 
results  of  the  integration  of  the  relevant  Variational'  equations  are 
compared  with  the  direct  integration  of  the  full  equations  of  motion. 

In  the  'hard'  forcing  cases  examples  are  given  showing  the  stability 
boundaries  in  the  excitation  amplitude/ frequency  plane. 

4.2  THE  TWO  MODE  SYSTEM 


Fig.  4.2.1 


By  approximating  the  experimental  structure  by  a two  mode  lumped 
mass  model  with  coordinates  as  shown  in  Figure  4,  the  equations  of 
motion  can  be  formulated  by  application  of  Lagrange's  procedure  to 
the  kinetic  and  potential  energy  functions.  The  components  of  'axial' 
motion  of  each  mass  (perpendicular  to  w^  and  are  kinematically 
related  to  the  end  displacements  w^  and  W2  by  the  assumption  of  a 


i 


deflection  curve  for  each  section,  similar  to  that  of  a cantilever  beam 
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under  static  end  load.  he  resulting  equations  of  motion  retaining  up  to 
quadrat ic  inert ial  nonl inearit ies  and  linear  parametric  terms  are  found 
to  bo: 


Imj  + + 2.25  jfilwj  - ‘ ‘"2  T7  ”2  ' "Tl 


1.5  m.,  + A,w, 

'1 


= -[mj  ^ m,(l  + 2.25  - 

..  ^2  *2  1.2  , *2a  0.3 

-m,l0.9  ^ + 0.45  (^2^2  ^ ^ '^l'^2 


3.0  ..  • • 

J-  + w^w^) 


1.5  m^  Wj  + + A^w, 

= '^2^s  ^ rf-  ^ '^l^P 


1.2,  ..  *2  ^ 

(WiW^  . wp  . 


1.2  .. 

" 17^2^ 


4.2.1 

where  A^  and  are  the  static  bending  stiffnesses  at  m^  and  m2  of  the 
cantilevers  of  length  and  I2  respectively,  which  make  up  the  model. 

4.3  NATURAL  FREQUENCIES  AND  NORMAL  MODES  - CONVERSION  TO  NORMAL 
COORDINATES 

The  equations  of  motion  4.2.1  may  be  written  in  matrix  form 
[A]w  + [C]w  = f(w,w,w,t)  4.3.1 

•St  ^ 


where  [A]  and  [C]  are  the  linear  symmetric  inertia  and  stiffness  matrices 
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and  f is  a column  vector  comprising  the  direct  forcing,  parametric 
and  quadratic  nonlinear  terms. 

The  homogeneous  part  of  4.3.1  yields  the  linear  natural 
frequencies  (eigenvalues)  and  normal  modes  (eigenvectors),  the  derivation 
of  which  may  be  found  in  any  standard  text  [8]. 

Applying  the  linear  transformation 

w = [R]p 

to  4.3.1,  ([R]  being  the  matrix  of  eigenvectors,  or  'modal'  matrix), 
we  have 


p + [R]-'[A]-'[C][R]p  = [R]''[A]''g(p,p,p,t)  4.3.2 

os#  ^ 

where  [R]” ^ [A]' ' [C] [R]  = [oj^]  is  the  diagonal  matrix  of  eigenvalues. 

Numerical  parameters  are  chosen  which  are  roughly  based  on  the 
experimental  system,  one  value  of  stiffness  being  adjusted  to  give  the 
required  internal  frequency  relationship. 

In  the  following  examples  the  values  of  the  various  parameters 
taken  are 

m^  = 0.3  kg,  m2  = 0.1  kg,  = 0.15  m,  = 0.1050  m 

= 2040.221  N/m,  A2  = 1.0  x 10^  N/m  4.3.3 

yielding  natural  frequencies  of  119.52  and  59.76  rads/sec. 

Viscous  damping  is  also  incorporated  and  again  rough  estimates 
based  upon  experimental  data  are  used  suggesting  values,  ^.d^j  = .012, 


"A 


j 

I 

I 

i 

1 
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Equations  4.3.2  are  equivalent  to  equations  3.1.1,  and  the 
coefficients  of  the  direct  forcing,  parametric  and  nonlinear  terms 
in  normal  coordinates  may  be  computed  through  4.2.1  and  4.3.3  via 
the  transformation  matrix  [R]"^[A]"^, 

4.4  'SOFT'  FORCING  RESULTS.  TWO  MODE  SYSTEM 

Figures  4.4.1  show  the  result  of  the  numerical  integration  of 
the  full  equations  4.2.1  and  the  corresponding  variational  equation 
3.4.1  when  excitation  is  at  the  lower  natural  frequency.  The  integration 
routine  used  in  all  cases  is  a fourth  order  variable  step  Runge  Kutta 
implemented  in  Fortran  on  a DEC  PDFS  minicomputer. 

The  points  are  plotted  to  the  scales  shown  such  that  the  actual 
values  are  given  by 

Pi  1...J  = Pi  » ,x  SC  1 + 0. 5 

‘ 1 plotted  ' 1 actual 

Pt  1 = p-  . , x SC2  - 0.5 

^2  plotted  ^2  actual 

where  SCI  and  SC2  are  given  on  the  figures. 

There  is  good  agreement  between  the  two  results  and  of  course  the 
integration  of  the  variational  equations  is  far  more  rapid.  Steady 
state  conditions  arc  achieved  and  so  it  is  possible  to  compute  the  response 
curves  for  this  resonance  case  by  solving  numerically  the  nonlinear 
algebraic  equations  which  result  from  setting  the  left  hand  sides  of 
3.4.1  to  zero.  Figures  4.4.2  show  some  results  of  this.  Here  to 


the  scales  shown  the  points  plotted  are  A(I)  x SC(I).  The  upper 
diagram  shows  the  results  obtained  for  the  complete  system  for  both 
cases  of  excitation  increasing  and  decreasing  from  the  lower  natural 
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Steady  state  response  curves  f2£  'v  - upper,  tuned  complete  system;  lower,  detuned  system.  (parametric  terms 
neglected. , 


frequency.  The  departure  from  the  linear  resonance  curve  for  mode  1 

is  apparent,  and  as  observed  experimentally  there  are  no  discontinuities 

in  the  curves.  It  is  also  noticeable  that  the  curves  are  not  quite 

symmetric  about  the  tuned  frequency,  and  it  has  been  found  that  removal 

of  the  parametric  (<  .)  terms  from  the  equations  of  motion  restore  the 

rj 

symmetry. 

The  lower  diagram  in  Figure  4.4.2  shows  the  effect  of  internal 
detuning.  The  parameters  of  the  system  are  adjusted  to  give 

= 1.01,  i.e.  becomes  59.171  rads/sec.  The  amount  of  detuning 
is  shown  on  the  diagram  and  here  the  parametric  terms  have  been  neglected. 
The  result  shows  a skewing  of  the  curves  about  the  lower  natural  frequency. 

Figures  4.4.3  show  the  result  of  numerical  integration  when  the 
excitation  is  at  the  higher  natural  frequency.  The  results  compare 
well  with  the  observed  experimental  behaviour  and  again  the  variational 
equations  3.4.2  model  the  complete  system  adequately.  Figures  4.4.4 
and  Figure  4.4.5  show  some  of  the  steady  state  response  curves  obtained 
for  this  case,  and  also  show  the  ordinary  linear  resonance  curve  of 
the  higher  mode.  In  the  tuned  system  symmetry  is  once  more  restored  by 
the  removal  of  the  parametric  terms  and  internal  detuning  (Fig.  4.4.5) 
again  causes  a substantial  skewing  of  the  response  diagram  about  the 
higher  natural  frequency.  The  outer  extremes  of  the  response  curves 
are  the  points  at  which  the  jumps  back  to  ordinary  linear  behaviour 
occur.  Unfortunately  the  resonance  region  could  not  be  approached 
from  the  'outside'  with  a view  to  finding  the  inner  boundaries  as  the 
solution  routine  used  will  not  converge  on  the  zero  solution  for  the 
lower  mode.  Analysis  [9]  shows  that  an  approximate  value  may  be  taken 
as  the  crossover  point  between  the  linear  and  nonlinear  response  curves 
for  the  higher  mode. 
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Fig.  4.4.4.  Steady  state  response  curves  'X-  - Upper,  complete  solution;  Lower,  solution  without  parametric  terms. 
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4.5  'HARD'  FORCING  RESULTS,  TWO  MODE  STRUCTURE 

Figure  4.5.1  shows  the  stability  diagram  obtained  for  the 
ji’’  combination  excitation  case,  + ^2  (see  equations  3.4.4  and 

3.6.11).  It  shows  the  effects  of  varying  the  damping,  removing  the 

■A 

D.C.  terms,  the  linear  forced  response  terms  (f^)  and  the  parametric 
terms  (i<^j.j)  separately  from  the  complete  analysis.  It  can  be  seen 
that  for  the  chosen  damping  values  (ed--  = 0.24)  excitation  amplitudes 
in  excess  of  30  mm  are  required  to  give  resonant  behaviour.  Physically 
this  is  a very  large  amplitude  and  does  not  compare  favourably  with  the 
experimental  model  which  exhibits  resonant  behaviour  above  3 mm 
excitation  amplitude.  However,  integration  of  the  full  equations  shows 
that  these  high  levels  are  indeed  required  to  initiate  resonant  behaviour. 

Hence  the  original  parameters  arbitrarily  chosen  for  the  mathematical 
model  are  presumably  not  sufficiently  precise  in  the  representation  of 
the  real  system. 

It  must  be  noted  that  neglect  of  either  the  parametric  (^^.j)  or 
the  linear  forced  response  terms  in  the  analyses  results  in  a much  wider 

i 

unstable  region  and  as  will  be  seen  shortly  gives  steady  state  results 

more  akin  to  the  experimental  observations.  ! 

The  upper  diagram  of  Fig.  4.5.2  shows  the  result  of  the  integration 
of  the  full  equations  when  excitation  is  within  the  unstable  or  interacting 
region,  the  lower  diagram  showing  an  expansion  of  the  first  two  seconds.  ] 

As  excitation  is  so  high  there  is  a large  component  of  forced  response 
on  each  mode,  but  it  is  clear  that  the  modal  frequencies  are  excited 
and  that  steady  state  conditions  are  not  achieved.  Figure  4.5.3  shows 
the  corresponding  complete  variational  simulation  at  these  excitation 
conditions  (equations  3.4.4)  and  the  various  components  of  the  response 
are  shown  separately  (note  different  scaling  factor  on  mode  2).  The 
different  transient  behaviour  is  caused  by  the  different  starting 
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conditions  for  each  solution,  but  the  'quasi'  steady-state  is  modelled 
reasonably  accurately  by  the  variational  analysis  except  for  an  error 
in  the  'beat'  period.  Figure  4.5.4  shows  that  neglect  of  the  D.C. 
terms  in  the  analysis  results  in  steady  state  conditions  and  an  unreasonably 
large  amplitude  for  the  lower  mode's  response. 

Figures  4.5.5  are  included  out  of  interest.  They  show  the  complete 
and  variational  simulation  for  this  case  when  the  parametric  terms 

are  neglected.  As  mentioned  above,  much  lower  excitation  is  required 
to  initiate  resonant  behaviour  and  the  results  compare  well  with 
experiment.  In  this  case  the  variational  simulation  models  the  full 
system  exactly. 

Figure  4.5.6  shows  the  stability  diagrams  obtained  at  twice  the 
higher  mode's  natural  frequency  (see  equations  3.4.7  and  3.6.13). 

Again  neglect  of  the  linear  forced  response  terms  or  the  parametric  terms 
results  in  a wider  region,  so  in  the  complete  solution  these  two  effects 
tend  to  negate  each  other.  Figure  4.5.7  shows  the  behaviour  of  the 
complete  system  when  excitation  is  within  the  unstable  region  (^^22  = .024). 
The  continuous  beating  between  the  modes  agrees  well  with  the  experimental 
results.  The  corresponding  variational  simulation  is  shown  in  Figure 
4.5.8.  Again  there  is  an  error  in  the  'beat'  period.  Here  neglect  of 
the  D.C.  terms  results  in  a 'beating'  exponential  growth  in  both  modes, 
and  so  the  D.C.  terms  are  essential  for  the  maintenance  of  the  reasonable 
accuracy  of  the  solution. 

Figure  4.5.9  shows  the  variation  of  the  respective  slowly  varying 
phases  with  time,  corresponding  to  Figure  4.5.8. 
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Figs.  4.5.5.  Full  and  variational  equations  with  parametric  (k^j)  terms 
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Fig.  4.6.2.  Arbitrarily  Displaced  System 

Figures  4.6.1  and  2 show  the  idealised  theoretical  model  chosen 
to  represent  the  four  mode  system.  Again  the  equations  of  motion  are 
derived  inclusive  of  the  quadratic  and  cubic  nonlinearities  via  formulation 
of  kinetic  and  potential  energy  functions  and  application  of  Lagrange's 
equations.  Static  deflection  curves  are  assumed  for  each  beam.  The 
resulting  equations  of  motion  in  the  four  generalised  coordinates 
Wp  V2>  ''’3r»  given  below  (Equations  4.6.1). 
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The  linear  homogeneous  equations  yield  the  natural  frequencies  and 
normal  modes  (Section  4.3),  and  these  have  been  found  numerically 
using  standard  matrix  solution  routines.  Taking  numerical  parameters 
based  on  experimental  measurement  of:- 


m^  = 144.35  gms 

= 1.3686  lO^gm/s 

^1  = 

135.5 

m^  = 143.10 

A^  = 3.67445  10® 

66.0 

m^  = 48.01 

Aj  = 9.28967  10** 

II 

97.0, 

the  resulting  frequencies  and  modal  matrix  is  shown  in  Figure  4.6.3. 
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It  can  be  seen  that  the  frequencies  obtained  are  not  exactly 
compatible  with  those  of  the  experimental  model.  As  expected  the 
higher  frequencies  are  overestimated,  and  so  the  desired  internal 
resonance  condition  (oj^  + oj,  is  not  exactly  fulfilled.  Indeed 

there  is  a better  internal  relationship  between  modes  2,  ^ and  4. 
However  inspection  of  the  nonlinear  terms  (transformed  to  normal 
coordinates)  shows  that  the  terms  required  to  give  resonant  behaviour 
(e.g.  ^9^4.  '”432’  simply  do  not  exist,  and  so  as  observed  in 

experiment  the  response  of  the  system  when  under  excitation  at  mode  4 
is  dominated  by  the  interaction  with  modes  1 and  3,  despite  the  large 
detuning. 

Figure  4.6.4  shows  the  result  of  the  direct  integration  of  the 
equations  of  motion  using  G.S.M.P.  simulation  language  on  an  IBM  360/50 
comjmtcr,  with  excitation  of  .07  mm  at  the  highest  natural  frequency. 
Here  'modal'  damping  based  on  logarithmic  decrement  tests  on  the 
exiierimenta 1 model  have  been  introduced  through  a Rayleigh  Dissipation 
Function  (stiffness/damping  ratio  of  30.10^).  The  result  clearly 
shows  interaction  between  the  excited  mode  and  the  antisymmetric  modes 
mucli  as  observed  in  experiment.  The  behaviour  of  the  other  mode  is 
interesting  as  it  exhibits  a D.C.  drift  with  ensuing  small  oscillations 
at  about  16  Hz.  Similar  behaviour  is  observed  experimentally  (Figure 
2.3.4).  This  phenomena  is  as  yet  not  understood  and  merits  further 
i nvestigation. 

Figure  4.6.5  shows  the  corresponding  result  obtained  from  the 
variational  equations  (equations  3.4.8).  There  is  good  correlation 
both  qualitatively  and  quantitatively.  It  was  found  to  be  impossible 
to  compute  steady  state  response  curves  for  this  case  as  a true  steady 
state  is  not  achieved.  The  slowly  varying  amplitudes  in  fact  exhibit 
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very  slight  variations  of  long  period  and  so  the  'stationary'  ((t)  = A = 0) 
ecpiations  have  no  solution.  A method  that  could  be  used  would  be  the 
repetetive  simulation  of  the  variational  eiiuations  followed  by  an 
estimation  of  the  average  steady  state  values  obtained. 

Figures  4.6.6  and  7 show  the  stability  diagram  and  a variational 
simulation  result  obtained  for  the  case  of  the  experimentally  observed 
combination  excitation  of  the  two  symmetric  modes,  leading  to  the  eventual 
participation  of  all  the  modes  of  the  system  (see  equations  .^.4.9,  3.6.11). 
In  fact  the  complete  system  was  found  to  have  an  infinite  stability 
boundary  in  this  case,  and  so  although  all  the  appropriate  terms 
required  by  the  equations  3.4.9  are  present,  this  is  a necessary  but 
not  sufficient  condition  for  resonant  behaviour  to  occur.  The  actual 


values  of  these  parameters  dictate  whether  or  not  the  system  will  be 
unstable.  The  figures  shown  were  obtained,  purely  for  the  purposes 
of  example,  by  neglect  of  the  linear  forced  response  terms  in  the 
analysis,  i.e.  the  system  is  reduced  to  one  having  simple  'parametric' 

terms  only.  Hence  once  again  the  theoretical  model  is  presumably 
not  sufficiently  precise  in  its  representation  of  the  experimental  system. 
Any  slight  modification  to  a system  parameter  will  result  in  changes 
to  the  modal  and  transformation  matrices,  which  in  turn  will  modify  the 
coefficients  of  all  the  nonlinear  terms,  and  hence  the  resulting  stability 
boundaries . 

Further  work  on  a simplified  version  of  the  four  mode  system  has 
been  carried  out  [9],  and  tlie  investigation  has  shown  that  inspection 
of  the  nonlinear  and  parametric  terms  present  in  the  governing  equations 
of  motion  indicates  which  of  all  the  possible  internal  and  external 
frequency  relationships  (section  3.5)  are  'safe',  simply  by  the  absence 
of  the  appropriate  terms  in  the  relevant  variational  equations.  Hence 
the  frequency  relationships  which  should  be  avoided  by  the  designer 
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Fig.  4.6.6.  Stability  Diagram,  Excitation  n£  near  (O2  + co^ 
(f.,  and  taken  arbitrarily  as  zero  to  give 
instability) 


Fig.  4.6.7.  Variational  simulation  when  within  unstable 
region  given  by  Fig.  4.6.6. 
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can  be  deduced. 
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However,  again  it  has  been  found,  particularly  in  the 
combination  situations,  that  although  the  relevant  terms  exist,  resonance 
might  still  not  occur,  the  values  of  the  coefficients  if  the  nonlinear 
terms  being  such  as  to  produce  an  infinite  stability  boundary. 


( 


E *■ 
E 

I'  ■ 


4 


-71- 


5. 


CONCLUSIONS 


Parametric  and  nonlinear  resonance  behaviour  in  structures  has 
been  illustrated  experimentally,  and  it  has  been  shown  that  quadratic 
inertia]  nonlinearities  are  of  primary  significance  in  determining 
the  response.  The  resonance  conditions  and  variational  equations 
listed  in  section  3 (and  Appendix  1)  can  be  adapted  and  applied  to 
any  multimode  structural  system  having  quadratic  inertial  nonlinearities 
and  linear  parametric  terms  subjected  to  harmonic  excitation.  The 
integrity  of  the  solution  technique  has  been  demonstrated  by  its 
application  to  the  two  and  four  mode  theoretical  models.  The  important 
factors  in  maintaining  the  accuracy  in  the  'hard'  forcing  cases  are  the 
inclusion  of  the  D.C.  terms  and  the  forced  response  components,  the 
latter  being  assumed  in  this  first  approximation  to  be  fixed  at  the 
linear  values.  The  inclusion  of  the  forced  response  terms  in  the 
generating  solution  give  rise  to  terms  of  nonlinear  origin  which 
accumulate  with  the  ordinary  linear  parametric  excitation  terms  in  the 
variational  equations.  The  combination  of  these  two  separate  effects 
may  either  reduce  or  increase  the  resonance  region  for  a given  structure 
depending  upon  the  values  of  certain  of  the  nonlinear  coupling  terms 
and  the  forced  response  of  al 1 the  modes  of  the  structure  (see  P.  j, 
equation  3.6.1). 

In  the  'soft'  forcing  cases  there  is  good  correlation  between 
the  theoretical  and  experimental  results  for  both  the  two  and  four 
mode  systems.  However,  in  some  of  the  'hard'  forcing  cases  (notably 


the  combination  case),  although  there  is  good  correlation  between  the 
full  equations  of  motion  and  the  corresponding  variational  equations. 
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comparison  with  the  experimental  results  is  relatively  poor.  This 
has  been  attributed  to  insut'fi  cient  ly  precise  representation  of  the 
experimental  systems.  However  in  these  cases,  particularly  in  those 
under  intense  excitation,  it  is  possible  that  the  cubic  nonlinearities 
may  also  be  important  in  determining  the  response,  as  the  resonance 
conditions  that  are  important  for  the  quadratics  are  also  important 
for  cubics  by  virtue  of  the  internal  resonance  conditions  (e.g. 

as  The  cubic  nonlinearities  can  easily 

be  incorporated  in  the  analysis,  but  this  will  naturally  lead  to  much 
more  complex  sets  of  variational  equations.  The  equations  of  motion 
of  most  space  structures  will  involve  the  quadratic  nonlinearities 
studied  here  and  in  general  they  will  be  the  most  important  nonlinear 
influence  on  the  response.  Higher  order  nonlinearities  will  however 
also  be  present  and  under  certain  circumstances  may  require  to  be  taken 
into  account.  Cubic  nonlinearities  will  introduce  further  possible 
external  and  internal  resonance  conditions  each  involving  up  to  three 
and  four  modes  respectively.  Consideration  of  even  higher  order 
nonlincarities  will  introduce  more  resonance  conditions  of  increasing 


complexity. 

The  experimental  work  indicated  the  form  of  the  relationship  between 
the  excitation  and  response  frequencies  in  multimode  systems.  The 
resonance  conditions  between  the  excitation  frequency  and  the  natural 
frequencies  (external)  or  amongst  the  natural  frequencies  themselves 
(internal)  need  not  be  met  precisely.  It  is  sufficient  that  they 
should  be  met  approximately,  the  response  frequencies  are  then  found 
to  be  related  to  the  exciting  frequency  in  such  a way  that  they  form  a 
linear  relationship  with  exact  integer  coefficients.  This  result  is 
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The  possibility  of  the  occurrence  of  'cascading'  type  resonant 
behaviour  has  also  been  shown.  It  may  be  concluded  that  for  more 
complicated  structures  having  a multiplicity  of  modes  and  particularly 
under  intense  excitation,  this  behaviour  may  be  exhibited  often. 

Indeed,  preliminary  experimental  investigations  into  thin-walled  box 
type  structures  under  intense  'parametric'  excitation  have  shown  that 
resonance  in  and  interaction  between  the  coupled  bending  and  cross- 
section  distortion  modes  is  quite  common,  and  in  general  the  resulting 
motion  involves  several  modes  and  is  noii-stationary . This,  along 
with  the  combined  direct  and  parametric  excitation  of  such  structures 
is  a promising  area  for  further  work.  The  other  area  of  immediate 
interest  for  further  research  lies  in  the  generalisation  of  the 
excitation  from  the  single  harmonic  considered  here  to  multifrequency 
excitation,  particularly  when  the  forced  and  parametric  excitations 
have  different  periods.  A single  mode  system  under  such  excitation  is 
considered  in  [10]. 
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APPENDIX  I i 

VARIATIONAL  EQUATIONS  FOR  THE  THREE  MODE  SYSTEMS  (Section  3.4.2)  i 

i 

j 

i 

1 , 1 5^^-  = + ^2  = ^3  (soft  forcing)  j 

j 

I 

+ F^cos4)^  - J [N^A^A^cos(<p^  - <p2  - <Pi)]  j 

-2A^n^  = F^sin(t)^  - j [-N^A^A^sin ((t)^  - <p2  - <P^)  - 


-2A2n2  = ' J [“^2^3^^'^  *■'*'3  ' ‘*’2^  ■ ^2^1^3®^'^*'‘^3  ‘ *^2  " ‘ ^^2‘^2'^22^2^ 

-2A,Q3$3  = (n^3  - a3pA3  - i [P3A2COS  (<{>3  - O)^)  - N3A^A2COS  ($3  - O)^  - 0^)] 
-2A3S  = - I [P3A2Sin(<f3  - (^2)  ^ N3A^A2sin ((1)3  - ct)^  - <))^)  - 2f23aj3d33A3] 


where 


^2  = 

^23 

^3  = 

*^32 

^ = 

‘^'"123  ■" 

"’l32^S^3  ■ 

£ 02 
123  3 

N2  = 

^'"213  " 

"'23p^lS  - 

‘^231^1 

^3  - 

■*-^312 

+ 11132^52^^2 

■ ^312^: 

0 q2 
132  2 

0 

213  3 

0 q2 
^32l“l 


1 , 2 + ^2  “ ^3  (soft  forcing) 


-2A^i2^$^  = (fi^  - wpAj  - I [P^A3COs(4)3  - + N^A2A3COS ((1)3  - <t>2  - 4>i)] 

-2A,fi,  = - J [-P^A  sin((t>3  - <))i)  - N^A2A3sin(<t)3  - 02  " " 2()jU>^dj ^A^ ] 
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- U)^JA,  + l■,C0S(j)2  - y [N2AjA^^cos(4)2  - (f^  - 0^)] 


-:A,i..  = I'-.sin^,  - y |-N^A^A_^.sin(0, 


- 2^2^21122^2] 


:A_i2_0,  = (.Qz  - CJ,)A,  - rr  [I'  A,  cos  {0  - 

.•)  .•)  .•)  J .•)  3 2 ‘ 3 1 3 


[J  + N^AjA2Cos(02  - 02  ■ 


■2A,!.2,  = 


- ^ [P,AjSin(03  - 0jJ  + N^^AjA2Sin(02  - 0,  - 0^3  - 


wlicrc 


1’=k:,,,  I\  = k_,  and  N , as  in  I.l. 
1 1 .•)  2 o 1 1,2,.') 


1 . 3 = 2Qj,  = Q,  (hard  forcing3 


-2A^i3|0j  = {n\  - w"^3A^  - j [P^AjCos20j  + N^A2A,^cos (0,  - 02  - 0^]  + 


.2AjQi 


- y [PjAjSin20j  - N^A^A3sin(0j  - 02  " 0i)  " 2fi^w^djjA^] 


-2A^Q,0,  = - uj^3A2  - j ' ‘]>2  ' '^2 


-ZA^Q-,  = - ^ [-N2AjA,sin(0,  ' ’^Z  ~ ' ^^2‘^2‘^22^2 ' 


-2A,fi30„  = (Q!  - w^3A^  - j [Ny\^A2Cos(02  - 02  - 0i)  I + ^3 


■2A,fi,  = 

J) 


- y [Ny\jA2sin{03  - (p^  - 0^ 


3 - 2fi,u),d,,A.] 
^ :>  3 33  .•)' 


where  N'  , are  as  in  I.l  and 

1 ) ^ ) J> 


‘’1  " '^11  ■ *^111^1 


P f 0^  - P f 0^  - P f 
^112  2“f  121  2“l  113  3 f 


-^3lS^l  " -"’111^1“^  " ^"'112  ^ "'121^^2^!%  " ^"’113  " ™131^^3S^f 


J 
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I • 4 = 29,^,  Uj  + ^2  = ^3  ('hard'  forcing) 


= ($2^j  - a)|)A^  - 2 [NjA^A^cosC^j  - ((.^  - (^jJ]  + 


2 [-NjA^A^sinCcJ)^^  - <^2  ' 'l>i)  ' 


“^2^2'*’2  ""  *-^2  ■ “P^2  " 2 ^2'^1^3^°'‘^ ^*^3  " ‘^*2  ' ^ 


■2\^n^  = 


- j [P2A2Sin2(()2  - N2A^A_sin ((l)^  - (p^  - 


2^3  = (^3  - - I [N3AiA2COs(03  - <1,2  - 0 ^ ) ] . V3 


-2'^3S  = - 2 tN3AjA2Sin(<^3  - (},2  - <^^)  - 2fi3U,3d33A3] 


where  = <22  - *'222^2  *^^2  ' '^1  *^^212^2  ^221^P  " ^3*^*'232^2 


*'223^f^  2m222f2S^f  *^™212  ™221^^1^2^f  ^"’223 


"'232^^3^2^f 


and  2 3 as  in  I.l 


1 .5  = 2f)3,  = ^3  (hard  forcing) 


■2N^i  = 


2 [-NjA2A3sin(<})3  - ' <*>1)  " 


l< 
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- 2 lN^A^A^cos((t)^  - 4>2  - 4>])]  + ^-, 


if 


- - [-N,A^A,siiU(})-  - <t>2  - 4)jJ  - 


-JA,Q,^„  = - ii)l).\„  - tr  [!'  A,cos2(}),  + N_A  A cos(i)),  - 4>9  - 4>,J]  + 

j .■)  .5  .•)  -•)  j)  2 3^  J 312  3 z 1 3 


■2AJ2,  = 

.■5  .■) 


\ lP,A,siii2cf),  + N-^A^A^siiUc}!,  - <{>2  " ~ ^^3‘^3‘^33'^-^ 


whore  N 


i ^ , are  as  in  I.l,  and 
1 , _ , 


p,  = K - + ai) 

J .•>3  .•>v*)J  .•)  1 


*^331  ^'l^f  ■ ^313^1^3  ■ ‘^332^2^f 
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Pj  - Kj2  - ^i22^2*'^2  ■ ^121^1^f  ‘ ^112^1^2  " ^123^3^f  ' 


”*'132^3^2  ■*■  * *-™121  "’ll2^^1^2^f  ^"'l23  "’l32^^3^2^f 


'’2  ■ '^21  ■ ^211^1*'^!  ■"  ■ ^212^2^f  ■ ^221^2^1  ' *'231^3^1 


4l3^'3^f  ■ ^"'211^1^1^f  ■ ‘■"’212  ■"  "’221^^2^1^f  ■ *■'"231  "'213^^3^1^f 


1.7  “ ^3  ^1  ^2  ~ ^3  forcing) 


- w^)A^  - y [N^A2AjCOs((i)2  - <t>2  ' <Pi)]  ^1 


^1^1  " ■ I [-N^A^A^sinCcf)^  - <p2  - 'P^)  " 2i2^U3^d^ ^A^ ] 


-2A2fi2^2  " ^^2  ■ ‘^^^2  ' I [P2A3COS  (c))^  + (fJj)  + N^A^A^cos  ((fj^  - 4>2  ' '}>^)]  + ^2 


■2A2^2 


[P^A^sinCcf.^  - P^)  - N2A^A3sin((t)3  - 4>2  ' <J>i)  - 2f22V22^2^ 


4S^3  = ^ - ‘^3^^3  - 2 [P3^2‘^°"^'^2  " ^3^  ^ ^3'^!^^°" ^<^3  ’ ^2  ' ‘^1^1  " ^ 


■2A3S  = 


2 [P3A2Sin((()2  ^ <(>3)  + N3A^A2sin ((^3  - ({)2  - (()j)  - 2fi3a)3d33A3] 


where  N are  as  in  I.l 

1 , Z f J 

^’2  " ^23  ■ ^213^1^3  ' ^231^1^f  ' ‘^233^3‘^^3  ^ ’ ^232^2^f 


■*^223^2^  ■"  ^"’232  ■"  '"223^^2^f^3  ^"'213  '"23p^l^3^f  ^"’233^3^3^f 


^3  " '^32  ■ *'312^1^2  ■ ^321^1^f  ‘ ‘*'322^2‘^^2  ' ^332^3^2 


■^323^3^f  ■"  ^"*332  '^323^  ^3^2^  f * ^’^312  * '"32p^l^2^f  ^‘^22^2^2^f' 


1.8 


- y [P^A^^cos((t)j  + (iJjJ  + N^AjAjCosCt})^^  - ^2  " ^ 


■-‘\^\  = 


j [PjA,sin((})j  + (})j)  - NjAjA^sin((t)^  - <t)^  - (p^)  - ^A^  ] 


-2A2fi24’2  = (^2  ■ *^^^2  ‘ 2 lN,A^A,cos  (<{),  - (J)^  - 0^]  + V, 


•2A,fi2  = 


- y [-N^A^A„sin(())2  - <P-^  - (p^)  - 2Q^iJ^d^^A^] 


■2AM^<p,  = (Q^  - - j iP^A^cosCcjjj  + (fij)  + NjA^A^cosCc})^  - (p^  - <))^)] 


" ■ y "■  NjA^A^sinCc]).^  - <t)2  - 4>j}  - 2Q^a)^d^^A^] 


where  N,  . , are  as  in  I.' 
1 , 2 , .•) 


and  P,  = K,,  - SL,  - £,  - £,,,£, 

1 Ij)  12^2j  132  2f  133  3 3 f 11.513 


-^3lfl^f  " ("'in  " "'l3l’^^fS  " ('"123  " "'132^  W3  " 


‘3  " ^31  ■ ^312^2^f  ■ ^321^2^1  " ^31 1^1  (^^1  ■"  " ^313^3^f 


+ (m,,,  + + (m,,_  + m,,  J f-fi.fl,,  + 2ni, , , f , 

.5.5 1 .5  1 .5.5 1 .513  3 1 f 312  321  2 1 r .5II  1 1 f 


